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ABSTRACT 

We examine the processes of formation and destruction of LiH molecules in the 
primordial gas at temperatures T < 5000 K in the framework of standard Friedmann 
cosmological models and we compute the optical depth of the Universe due to elastic 
Thomson scattering of cosmic background photons on LiH. With the help of a simple 
model of evolution of a spherical density perturbation we follow the linear growth 
and the subsequent collapse of protogalactic clouds of various masses, evaluating con- 
sistently their physical characteristics and chemical composition. We determine the 
expected level of anisotropy of the cosmic background radiation due to the presence 
of LiH molecules in high-redshift protogalactic clouds. We conclude that LiH spectral 
features generated in the primordial gas are hardly detectable with current techniques 
and instrumentation. 

Key words: early Universe - molecular processes. 



1 INTRODUCTION 



of secondary anisotropies in the CBR temperature at a fre- 
quency f is simply 



It has been suggested that a finite amount of molecules 
such as Ha, H+, H+, HD, HD+, H2D+, HeH+, LiH and 
_ ■ LiH^ in the primordial gas can induce in the cosmic back- 
, ground radiation (CBR) both spectral distortions and spatial 
» I anisotropies (see Dubrovich 1994 and Maoli 1994 for recent 
' reviews). For example, if the CBR spectrum had some in- 
trinsic distortion in the Wien region originated at high red- 
shift (2 ^ 100) by whatsoever process of early energy re- 
lease, then molecules act as heat pumps, absorbing energy 
in their rotovibrational levels and re-emitting this energy at 
lower frequencies via rotational transitions (Dubrovich 1977, 
1983). Clearly, such a process can occur only if the CBR 
spectrum deviates from a pure Planck profile. Dubrovich & 
Lipovka (1995) have performed detailed calculations of this 
effect for HaD"*" molecules: the distortions appear like narrow 
absorption and emission features in the CBR spectrum. 

Spatial anisotropies in the CBR can be produced by 
Thomson scattering of photons on molecules (or electrons) 
located in protoclouds moving with a peculiar radial ve- 
locity component Vp^c (Sunyaev & Zel'dovich 1972, 1980; 
Dubrovich 1977, 1983; MaoU 1994; MaoU, Melchiorri & 
Tosti 1994). If the abundance of molecules in the primor- 
dial gas is sufficiently high, the cumulative effect of elastic 
Thomson scattering can even create a "curtain" that blurs 
primordial CBR anisotropies (Dubrovich 1993). The level 
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where Ti, is the optical depth at frequency f through the 
cloud. The advantage offered by molecules is that they pos- 
sess much larger photon scattering cross sections (by more 
than 10 orders of magnitudes) than free electrons. The 
derivation of the resulting AT/T requires the solution of 
the equation of radiative transfer in the expanding Universe, 
described in detail in Appendix A. 

In addition, the presence of molecules in protogalactic 
clouds may induce small ffuctuations in the CBR tempera- 
ture when the population of rotovibrational levels deviates 
from the equilibrium (Boltzmann) distribution correspond- 
ing to the actual radiation temperature T^. Since the tem- 
perature of the gas Tg and the radiation temperature are dif- 
ferent for z ^ 10'', the populations of molecular levels may 
differ from a Boltzmann distribution with temperature Tr in 
clouds of sufficiently high density where coUisional excita- 
tion processes become significant. The level of of anisotropy 
produced in this case is proportional to the difference be- 
tween Tr and the excitation temperature of the molecular 
transition. 

The efficiency of the coupling between CBR photons 
and the primordial gas clearly depends on the chemical 
composition of the gas, and therefore a description of the 
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Figure 1. The radiative association rate i?ass of LiH as func- 
tion of the gas temperature Tg. The dashed line shows our fit to 
the values computed by Dalgarno et al. (1996) and Gianturco & 
Gori Giorgi (1996a) (squares); the dot long dashed and dot short 
dashed line our fit to the numerical results by Gianturco & Gori 
Giorgi (1996b) (solid and empty circles) for the B^Il X^S"*" 
and A^T,'^ X-'E+ processes respectively. The relative popu- 
lations of the ^5 and levels of Li have been assumed to be 
at thermal equilibrium with temperature Tr = Tg. The solid line 
shows the total rate. 



chemistry of the early Universe is necessary. A large number 
of papers has been devoted to this problem (e.g. Dalgarno 
& Lcpp 1987; Palla 1988; Puy et al. 1993; Palla, Galh & 
Silk 1995). 

Of all primordial molecules, lithium hydride (LiH) is of 
considerable interest since the high dipole moment (5.9 de- 
byes) of this molecule makes its rotational and rotovibra- 
tional transitions particularly strong (see e.g. Zemke & 
Stwalley 1980, Bellini et al. 1994, Gianturco et al. 1996). 
For this reason, the possibility of detecting LiH features in 
the spectrum of the GBR and in high redshift protogalac- 
tic clouds has raised a considerable interest both theoreti- 
cally (Dubrovich 1994, MaoU 1994, Maoh et al. 1994, MaoU 
et al. 1996) and observationally (de Bernardis et al. 1993). 
Specific observational programs are in progress (Signore et 
al. 1994). 

It is important to remind that the abundance of LiH in 
the pre-galactic gas is proportional to the primordial abun- 
dance of Li (hereafter [Li/H]p), which is still highly contro- 
versial. It was originally shown by Spite & Spite (1982) that 
the Li abundance on the surface of very metal-poor Pop 
II stars ([Li/H]pop ii — 10^^") is basically independent on 
metallicity (the so-called Spite plateau). The Pop II value 
compares well with the abundances of the other light ele- 
ments predicted by standard big-bang nucleosynthesis and 
might therefore represent the true primordial Li: [Li/H]p ~ 
[Li/H]pop II ~ 2 X 10~^°. However, evidence is accumulating 
for a higher primordial value, of the order of [Li/H]p ~ 10~^. 
Prom observations made with the Keck telescope of a sam- 



ple of sub-giants in the globular cluster M92, Deliyannis, 
Boesgaard & King (1995) have recently found evidence of 
significant dispersion in the measured values of [Li/H]pop ii 
values, strongly suggesting that some amount of depletion 
has occurred among Spite's plateau stars. If this were the 
case, then the Li abundance of Pop II stars could be the 
result of depletion from a higher initial primordial value, a 
result consistent with the rotating stellar models by Pinson- 
neault, Deliyannis & Demarque (1992). As a possible way 
to distinguish between the high- [Li] p and the low- [Li] p sce- 
narios, Steigman (1996) has recently proposed to determine 
the ratio of [Li/K] from absorption features along differ- 
ent lines-of-sight in the Galaxy and the Large Magellanic 
Gloud. At the present stage, available data along the line- 
of-sight of SN1987A favor a primordial abundance of Li 
lower by at least a factor ~ 2 than the present Pop I value 
([Li/H]pop I ~ 1.6 X 10-^). 

Since the main quantities determined in this study 
(abundances, optical depths and anisotropies) depend on 
the assumed initial Li abundance, we should expect an over- 
all increase of their values by an order of magnitude, if the 
"high" value of [Li]p were to be confirmed by future obser- 
vations. 



2 FORMATION AND DESTRUCTION OF 
LITHIUM HYDRIDE 

The LiH molecule has been the subject of a large number 
of theoretical investigations aimed at determining the con- 
ditions for its formation in various astrophysical contexts. 
Kirby & Dalgarno (1978) (see also Bochkarev & Kherson- 
kii 1985) considered the photodissociation of LiH following 
absorption from the v — vibrational level of the ground 
X^'S'^ state into the vibrational continuum of the A^S"*" 
or B^H excited states. The cross sections and rates for the 
inverse radiative association processes 

Li(^P) -I- H ^ LiH(yl^E+) ^ LiH(X^E+) -|- hu, 

uCP) + H ^ LiH(B^n) ^ LiH(X^E+) + hu, 

have been recently computed by Gianturco & Gori 
Giorgi (1996b). Since the Li(^P) level lies about 1.85 eV 
above the Li(^S') level, the formation of LiH molecules via 
electronically excited levels is important only when the tem- 
perature is greater than some thousand degrees K. At lower 
temperatures, processes involving transitions between vibra- 
tional levels of the ground X^E"*" state and its continuum 
become the dominant mechanisms of radiative association 
and dissociation of LiH. Such physical conditions are real- 
ized in the interior of dense clouds, and in the primordial 
pre-galactic gas at epoch z ^ 500. 

For the direct radiative association of Li(^S) and H, 

Li(^S') -h H -> UIl{X^T,+) + hv, 

Lepp & ShuU (1984), using semi-classical arguments, esti- 
mated a rate coefficient Rasa — 10~^^ cm^ s~^, a value 
adopted in several studies of the chemistry of the early Uni- 
verse (e.g. Lepp & ShuU 1984; Puy et al. 1993; Palla et 
al. 1995). Khersonskii & Lipovka (1993) performed a fully 
quantum mechanical calculation analysis of this reaction. 
Their results show a rapid increase of the cross sections for 
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Table 1. List of considered reactions 





reaction 




rate (cm^ s ^) 


notes 


reference 


1) 


Li+ + e — > Li + /ii/ 




1.036 X 10-"[i/Tg/107.7x 
(14- ^Tg/107.7)«-e6i2x 












Cl 4- -kIT^IX 177 X io71i-3388i-i 


miantal palpnlatinn 


yt%) 




lA^-hv ^ Li+ + e 




detailed balance applied to (1) 




3) 


Li(2 5) + H -» LiH(XiE+) + hv 


10-" 


semiclassical estimate 


(b) 








1.0 X 10-i6Tg-o «i 


quantal calculation 


(c) 








(5.6 X 10^'^Tg-"'^^ 4- 7.2 x lO'^-'T^-'^'^)-'^ 


quantal calculation (adopted) 


{d),(c) 


4) 


LiH(X^S+) hv ^ Li(25) - 


f H 




detailed balance applied to (3) 




5) 


Li(2p) + H -> LiH(XiS+) 4 


- hv 


2.6 x 10-i'*TO i2cxp(-rg/6500) 


AiS+ XiS+, quantal calc. 


(f) 


6) 


LiH(XiE+) — ► Li(2p) + H 4 


- hv 




detailed balance applied to (5) 




7) 


Li(2p) + H LiH(XiS+) 4 


- hv 


1.9 X 10-"Tg-°-34 


fiin -* XiS+, quantal calc. 


(f) 


8) 


LiH(XiS+) ->■ Li(2p) + H 4 


- hv 


6.3 X 10-''Tg~^^^ - 7.6 X 10-^4- 

9 V 1 n — IOT"*^/^ 

Z.U A ±U J- g 


detailed balance applied to (7) 






Li+ 4- H- -> Li 4- H 




■ 

expemiientai cross section 






Li 4- e — > Li" 4- hv 




5 7 X 10-l'''T°-5^ pxnC— T„/17200'l 


experimental cross section 




111 


\A~ 4- /ii^ — > Li 4- e 




detailed balance applied to (10) 




12") 


Li 4- H+ — ► Li+ 4- H 




2 'i X in-*''7-'7.9p„_/_T-' /-loin") 

^ .0 A J-tJ S i^V S / J 








Li + H+ ->■ Li+ 4- H 4- hv 




1 7 in — ISt"— 0.051 rji /oQonnn^ 
i.( X lU Jg 6Xp^^ — Jg/z<5zUUUj 


quantal calculation 


(J) 


14j 


Li~ 4- H+ ^ Li 4- H 






same as (9) 


(g) 


15) 


Li+ + H -> LiH+ + hv 




A A in — 20^ — 0.9 „..„/ ^ /'7nnn\ 
1.4 X lU ig^ exp(— ./g/7UU0j 


quantal calculation 


(a).(e) 


16) 


LiH+ 4- fti^ ^ Li+ 4- H 






detailed balance applied to (15) 




17) 


LiH+ + c ^ Li + H 




47 

3.8 X 10 Tg ■ 


estimate 


(k) 


18) 


LiH+ 4- H ^ Li 4- H+ 




9.0 X 10-1" exp(— 66400/Tg) 


estimate 


(k) 


19) 


LiH+ 4- H -» Li+ 4- H2 




3.0 X 10-1° 


estimate 


(k) 


20) 


LiH+ 4- H ^ LiH 4- H+ 




1.0 X 10-" exp(-67900/Tg) 


estimate 


(k) 


21) 


Li 4- H- ^ LiH 4- e 




4.0 X 10-10 


estimate 


(k) 


22) 


Li- + H ^ LiH 4- c 




4.0 X 10-1° 


estimate 


(k) 


23) 


LiH + H+ ^ LiH+ 4- H 




1.0 X 10-9 


estimate 


(k) 


24) 


LiH + H ^ Li 4- H2 




2.0 X 10-11 


estimate 


(k) 


25) 


LiH + H+ ^ LiH+ 4- H 




1.0 X lO-'J 


estimate 


(k) 


26) 


LiH + H+ ^ Li+ + H2 




1.0 X 10-" 


estimate 


(k) 


27) 


LiH 4- H+ ^ Li+ 4- H+ 




1.0 X 10-'' 


estimate 


(k) 



(a) Verner & Ferland (1996); (b) Lepp & Sliull (1984); (c) Khersonskii & Lipovka (1993); (d) Dalgarno, Kirby & Stancil (1996); 
(e) Gianturco & Gori Giorgi (1996a); (f) Gianturco & Gori Giorgi (1996b); (g) Peart & Hayton (1994); (h) Ramsbottom, Bell & 
Berrington (1994); (i) Kimura, Dutta & Shimakura (1994); (j) Stancil & Zygelman (1996) (k) Stancil, Lepp & Dalgarno (1996). 



radiative association with the vibrational quantum number 
v. the molecule is formed preferencially in vibrational levels 
with V =15-18. The resulting rate coefHcient for the direct 
radiative association of LiH is in good agreement with the 
semi-classical estimate of Lepp & ShuU (1984) (see Table 1). 

More recently, the LiH chemistry in the primordial gas 
has been reanalyzed in detail by Dalgarno, Kirby & Stan- 
cil (1996) and by Stancil, Lepp & Dalgarno (1996, hereafter 
SLD) . The accurate quantal calculation of the radiative asso- 
ciation rate of LiH performed by Dalgarno et al. (1996) gives 
a value of i?ass smaller than the one computed by Lepp & 
ShuU (1984) and Khersonskii & Lipovka (1993) by two-three 
orders of magnitude. The reason for such a large diflference 
is not clear, but it may be related to the approximations 
adopted by the latter authors. The quantal calcualation of 
J?ass for reaction (4) recently performed by Gianturco & Gori 
Giorgi (1996a) has confirmed SLD's results to a very high 
degree of precision. This is even more remarkable, since the 
two calculations arc based on independent sets of potential 
curves and wavefunctions. 

The total rate of radiative association of LiH is repre- 
sented graphically in Fig. 1 as function of the temperature 



Tg for a gas where the relative populations of the S and 
levels of Li are assumed to be in thermal equilibrium at 
temperature Tg. It is clear from the Figure that radiative 
association via excited electronic levels is the dominant for- 
mation process of LiH when the gas temperature is larger 
than ^ 2000 K. Otherwise, LiH is formed only by direct 
radiative association. 



3 CHEMICAL EVOLUTION OF THE 
PRIMORDIAL GAS 

We consider now the chemical evolution of the pregalactic 
gas in the framework of Fricdmann cosmological models. 
Our standard model is characterized by a Hubble constant 
Hq = 67 km Mpc-^ (van den Bergh 1989), a closure pa- 
rameter f2o = 1, a baryon-to-photon ratio r?io = 4.5 (Galli et 
al. 1995), and a present radiation temperature To = 2.726 K 
(Mather et al. 1993). The initial fractional abundance of 
Li nuclei for this model is 2.1 x IQ-i" (Smith, Kawano & 
Malaney 1993). The integration of the chemical network 
starts at z = 2500, with H and Li fully ionized and He fully 
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Table 2. Final fractional abundances 





Vio = 0.7 








T]io = 4.5 








mo = 10 








z = 10 




z = 




z = 10 




z = 




z = W 




z = 




H+, e 


4.6 X 10" 


3 


4.6 X 10" 


3 


7.2 X IQ- 


4 


7.1 X 10" 


-4 


3.2 X 10- 


4 


3.2 X IQ- 


4 


H- 


1.9 X 10" 


11 


1.9 X IQ- 


13 


1.5 X IQ- 


12 


1.5 X 10" 


-14 


4.9 X 10- 


13 


4.9 X 10- 


16 


Li 


8.2 X 10" 


10 


8.5 X 10- 


10 


1.6 X IQ- 


10 


1.7 X 10" 


-10 


8.9 X 10- 


10 


9.3 X IQ- 


10 


Li+ 


2.8 X 10- 


10 


2.5 X 10" 


10 


4.5 X IQ- 


11 


3.8 X 10- 


-11 


2.1 X 10- 


10 


1.7 X 10- 


10 


Li- 


5.0 X 10- 


20 


4.2 X 10" 


22 


4.3 X IQ- 


21 


3.9 X 10- 


23 


1.6 X 10- 


20 


1.5 X IQ- 


22 


LiH 


1.5 X 10- 


18 


4.2 X 10- 


19 


2.2 X IQ- 


19 


9.8 X 10- 


20 


1.2 X 10- 


18 


4.8 X 10- 


19 


LiH+ 


8.1 X 10" 


17 


1.2 X 10" 


16 


8.5 X IQ- 


18 


2.2 X 10- 


17 


2.9 X 10- 


17 


1.0 X 10- 


16 




log (1-l-z) 




Figure 2. Fractional abundance (by number) of atoms, ions and 
molecules as function of the rcdshift for the standard model. The 
effect of the radiative association of LiH in excited electronic levels 
is shown by the dashed line. 



Figure 3. Processing rates of the most important reactions for 

the formation of LiH. The curves arc labelled with the corre- 
sponding reaction numbers listed in Table 1 



recombined. We have also considered the possibihty of a high 
[Li/H]p, as discussed in the Introduction. If [Li/H]p is close 
to 10~^, standard big-bang nucleosynthesis allows either a 
"low" or a "high" value of 7710 (7)10 — 0.7 or ~ 10, respec- 
tively). The primordial Li abundance predicted by models of 
inhomogeneous Big Bang nucleosynthesis can even be con- 
siderably larger than the values adopted here, but since its 
value depends sensitively on a number of parameters and 
assumptions (see e.g. Mathews et al. 1990), for simplicity 
we consider models based only on standard Big Bang nucle- 
osynthesis results. The initial fractional abundances of H, 
He and Li are 0.936, 0.064 and 1.1 x 10"'' for the rj = 0.7 
model, 0.926, 0.074 and 2.1 x IQ-^" for the rj = 4.5 model, 
and 0.923, 0.077 and 1.1 x 10"'' for the 77 = 10 model. 

Our chemical network includes about one hundred re- 
actions among 21 chemical species. Reactions involving Li- 
boaring species and the adopted rates arc listed in Table 1. 
For the recombination rate (reaction 1) we have adopted the 
accurate fitting formula given by Verner & Ferland (1996), 
based on the photoionization cross sections computed by Li 



by Peach, Saraph & Seaton (1988) and Gould (1978). The re- 
combination rate coefficient of Li shown in Table 1 is in good 
agreement with the results by Caves & Dalgarno (1978). 

The evolution of the gas temperature Tg is governed by 
the equation (see e.g. Galli 1990; Puy et al. 1993; Palla et 
al. 1995) 

dr Ft 2 

'W = ~'^'^^R + Sfc^^t^^ " ^)Compton + (P - A)^ol]. (2) 

The first term represents the adiabatic cooling associated to 
the expansion of the Universe {R is the scale factor of the 
Universe). The other two terms represent respectively the 
net transfer of energy from the CBR to gas (per unit time 
and unit volume) via Compton scattering of CBR photons 
on electrons, 

rr A^.. _ 4fccTTaT,^(T, - Tg) 

11 — AjCompton — Tie, (o) 

TTleC 

and via excitation and de-excitation of molecular transi- 
tions. 
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(r - A)moi = ^(rifCij - njCji)hvij, (4) 

where Cij and Cji are the coUisional excitation and de- 
excitation coefficients and rii are the level populations. The 
adopted expressions for the radiative and collisional coeffi- 
cients of LiH as well as the procedure followed to determine 
the level populations are summarized in Appendix B. For 
the molecular heating and cooling of the gas we have con- 
sidered the contributions of H2, HD and LiH. 

When the rate of collisional de-excitation of the roto- 
vibrational molecular levels is faster (slower) than their ra- 
diative decay, the energy transfer function (F — A)moi can 
become an effective heating (cooling) source for the gas (cf. 
Khersonskii 1986, Puy et al. 1993). In general, the contribu- 
tion of molecules to the heating/cooling of the gas becomes 
important only during the phase of collapse of protogalactic 
objects (see Sect. 4). 

The chemical/thermal network is completed by the 
equation for the redshift 

^ = 1 (5) 

dz Hq(1 + z)'W^ + ^oz' 

where Hq is the Hubble constant. The density n{z) of 
baryons at redshift z is 

n(«) = f2bficr(l -I- a)^, (6) 

where ricr is the critical density. The abundance of atoms 
and molecules is measured with respect to the total number 
of baryons: /(H) = n(H)/n and so on. 



/ \ 
/ \ 
/ \ 
/ \ 
/ \ 




1 10 100 1000 

frequency (GHz) 



Figure 4. Optical depth of the Universe generated by Thom- 
son scattering of CBR photons on LiH molecules for the standard 
model {solid curve) and for the high primordial Li cases {dot short 
dashed curve, rjio = 0.7; dot long dashed curve, 7710 = 10). For 
the standard model, the contribution of rotational and rotovibra- 
tional transitions is shown by the dashed curves. 



3.1 Results 

The evolution of the fractional abundances of H, H"*", H~, 
Li, Li^, Li~, LiH^ and LiH for the standard model is shown 
in Fig. 2 as function of the redshift z. The abundances at 
z = 10 and 2; = for all the models considered are listed in 
Table 2. 

In all models, the final abundance of LiH is extremely 
small, of the order of 10"^''-10"^®. The maximum final 
abundance (~ 5 x 10~^^) is obtained for the 7710 = 10 
model. Fig. 3 shows the relative contribution of individual 
chemical reactions to the formation of LiH in the early Uni- 
verse. Radiative association from Li('^P) is significant only 
for z > 650, with the process B^H X^T,'^ being dom- 
inant. Otherwise, LiH is formed mostly by radiative asso- 
ciation from Li(^S') (reaction 3). At « ~ 80 the reactions 
of associative detachment (21) and, most importantly, (22) 
also contribute to the formation of LiH. Notice however that 
although the value of the rate coefficients for the latter re- 
actions is rather uncertain (see SLD), the asymptotic value 
of the abundance of LiH is largely determined by reaction 
(3), whoso rate is now well established (Dalgarno, Kirby & 
Standi 1996, Gianturco & Gori Giorgi 1996a). 

The evolution of Li~ for z < 100 closely mimicks that 
of H~ , both ions being removed at low temperature by re- 
actions of mutual neutralization with H"*" at a rate propor- 
tional to Tg . Since Tg (X (1 + z)'^ in this range of rcsliift, 
the abundances of H~ and Li^ decrease like (l+z), as shown 
in Fig. 2. Our computed abundances of H"*", H^, Li and 
LiH"*" at z = 10 are in excellent agreement with the results 



obtained by SLD (within 10-20%), our abundance of Li"*" be- 
ing a factor 1.8 smaller. Our values for Li~ and LiH at a = 10 
are smaller by large factors (~ 900 and ~ 300, respectively) 
than those published by SLD, but the two computations give 
essentially the same results when a sign mistake in the rate 
adopted by SLD for reaction (10) is corrected. 

An interesting result is that LiH"*" is more abundant 
than LiH for 2 ^ 18, by a factor ~ 200 for the standard 
model. This is in agreement with the early prediction by 
Dalgarno & Lepp (1987). 



3.2 Optical depth 

Having obtained the LiH abundance as function of the red- 
shift, we are now able to compute, with the help of the for- 
mulae listed in Appendix A, the total optical depth of the 
Universe due to elastic Thomson scattering of CBR photons 
on LiH molecules (see eq. [A9]). The result is shown in Fig. 4 
as function of the frequency of observation z/. The maximum 
optical depth for the standard model is rf^^ii = 6 x 10~^^, 
at p"^^^ — 30 GHz. Pure rotational transitions of the v = 
level dominate at this frequency. The general behaviour of 
rLiH(i^) shown in Fig. 4 is in agreement with the one obtained 
by Maoli et al. (1994). The attenuation of CBR anisotropics 
being of order tlih, little room is left for the "molecular 
blurring" suggested by Dubrovich (1993). Even for the high 
[Li/H]p model with rjio = 10, where a larger optical depth 
is obtained, the induced attenuation is still insignificant. 
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Table 3. Physical (]naiititios for a lO'-" collapsing jji-otogalaxy 



t/tfr 


z 


71 




g 


R 


e 




V 


J \ -^J 




f(UH) 








( cm~^^ 




(K) 


(knc) 


V / 


(km s"''") 


(km s^^) 










0.0 


10.00 


1.38 X 10- 


-3 


9.6 


38.8 


28.2 


0.0 


90 


4.5 X 10-*' 


5.3 X 10-'^ 


2.2 X 10- 


19 


0.1 


7.45 


1.41 X 10" 


-3 


9.8 


38.5 


22.8 


3.8 


103 


4.5 X 10-^ 


5.4 X 10-'^ 


2.4 X 10" 


■19 




D.Uo 


1 /in \y 1 n~ 


-3 




O ( .o 


1 n 

iy.o 


■7 

( .0 


1 1 Q 

iio 


A a \y ^ n — 6 
4.D X lU 


0.4 X lU 


^.o X iU 


■19 
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4 A SIMPLE MODEL OF A COLLAPSING 
PRIMORDIAL CLOUD 

So far our attention has been directed exclusively to the 
global effects of LiH molecules in the homogeneous primor- 
dial gas. One of the main limitations to the occurence of de- 
tectable features is the low molecular abundance obtained 
in such a context. It is therefore of interest to follow the 
chemical (and thermal) evolution of the first density inhomo- 
geneities that emerged from the expanding primordial gas. 
Those initial density perturbations whose mass exceeded 
the Jeans mass were unstable to gravitational collapse and 
evolved to form galaxies and galaxy clusters. The increase in 
density during the first (quasi adiabatic) phase of collapse 
can provide suitable conditions for the build-up of consider- 
ably larger molecular fraction than in the surrounding ex- 
panding Universe. In addition, velocity gradients in a cloud 
may contribute to enhance the effective optical depth of lo- 
calized regions in the Universe (Maoli et al. 1994, 1996). 

Let us consider then the collapse of a purely baryonic 
primordial cloud in order to determine the abundance of the 
relevant molecular species as function of time and to predict 
the level of anisotropies induced by molecular transitions 
during the collapse of the cloud. We follow the approach of 
Lahav (1986) considering a spherical homogeneous cloud of 
meiss M, uniform density n and radius R contracting (or 
expanding) in a homologous fashion (v(r) oc r). We assume 
that the radius of the cloud reaches its maximum value at the 
turn-around epoch Ztz, when the density and temperature of 
the cloud are related to the corresponding quantities in the 
ambient gas by 



Tg(«ta). 



(8) 



The radius at turn around i?ta is fixed by the mass and 
density of the cloud. In adiabatic collapse the temperature 
of the gas in the cloud, Tg, increases like R~^. 

By virtue of the hypothesis of homologous collapse, the 
primordial cloud may be treated in the same way as the 
expanding Universe in the previous Sections. For example, 
eq. (2) for the gas temperature still holds if one replaces 
the scale factor of the Universe with the cloud's radius. The 
optical depth of the cloud by elastic Thomson scattering of 



CBR photons on LiH molecules is given again by eq. (A9), 
but the cosmological velocity gradient is now replaced by 
the velocity gradient due to gravitational collapse. A more 
sophisticated calculation of the optical depth, taking into ac- 
count projection effects, can be found in Maoli et al. (1996). 
Defining a non-dimensional radius ^ and temperature 6 

by 



R — ^Rti 



(9) 



the equations for the radius and for the temperature can be 
written as (Lahav 1986) 



d^$ _ 5fcTt 



-3 GM^^_2 



(10) 



dt^ UniHR^, 
and 

d^ ^ 2 (F - A)compton + (r - A)^ol 2 

dt 3 nfcTta ' 

It is convenient to express time in units of the free-fall time 
for the cloud, 



(11) 



3 \ 1/2 



(12) 



Once the mass Al and the turn-around rcdshift zta are fixed, 
the evolution of the model cloud is determined by the in- 
tegration of eq. (10) and (11). A relation between M and 
2ta can be obtained by prescribing the power spectrum of 
the density perturbations (see e.g. Gunn & Gott 1972, La- 
hav 1986). For simplicity, we prefer to use a reference value 
2ta = 10 for all masses. 

We consider here a low-mass cloud with M = 5 x 
10^ Mq, slightly above the Jeans mass at Zta. = 10 (Mj — 
3 X 10^ Mq), and a galactic-size cloud with M = 10^° Mq. 
Fig. 5 shows the time evolution of the cloud's radius, den- 
sity and temperature (all normalized to ther initial turn- 
around values) for the M = 5 x 10^ M© case. The evolu- 
tion is followed for about three free-fall times for illustrative 
purposes only; in a more realistic model, relaxing the as- 
sumption of homologous contraction/expansion, the actual 
evolution of the cloud is expected to deviate from the one 
shown here owing to the shocking of different mass shells 
(see e.g. Hairnan, Tlioul & Loeb 1996). Because of molecular 
cooling, the model cloud performs a series of oscillations of 
decreasing amplitude and period, in striking contrast with 
the corresponding evolution in the adiabatic case (dashed 
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Figure 5. Time evolution of radius, density and temperature for 
a primordial cloud of mass M = 5 x 10^ Mq. All quantities are 
normalized to their initial turn around values. Time is in units of 
the free-fall time. Dashed lines: adiabatic evolution; solid lines: 
non-adiabatic evolution. 




<; 

I 



-LiH xlO^ 



H — [ — \ — I — ! — h 




Figure 6. Upper panel: fractional abundances of the relevant 
atomic and molecular species as function of time for the same pri- 
mordial cloud as in Fig. 5. Lower panel: absolute value of the net 
energy transfer due to molecular and Compton heating/cooling 
processes. 



lines). Fig. 6 shows the fractional abundance of the relevant 
atomic and molecular species and the net energy exchange 
between gas and radiation due to molecular and Compton 
processes, for the same case as in Fig. 5. Both processes act 
as a net heating for the cloud for t < 0.53 tg, when Tr > Tg, 
and as a net cooling afterwards. It is important to notice the 
relevant role played by HD molecules, in addition to H2, for 
cooling the gas during the early phases of collapse. As shown 
in Fig. 6, heating/cooling by LiH molecules is negligible at 
all times. 

The relevant physical characteristics during the collapse 
of the M = Mo cloud are summarized in Table 3. The 
listed quantities arc: rcdshift 2, cloud density n, temperature 
Tg, radius R, angular diameter 9, collapse velocity iicoii of the 
outermost shell, peculiar velocity Upcc, and fractional abun- 
dances of the most relevant molecular species. The angular 
diameter of a spherical cloud of proper radius R at redshift 
z is computed from the formula (see e.g. Peebles 1971) 



e{R,z) = 



HoR 



noz + {no - 2)(vi + no2- 1)' 



(13) 



For the peculiar velocity we follow Maoli et al. (1994) as- 
suming the simple law of evolution of velocity dispersion in 
the linear regime (see e.g. Peebles 1980) 



= (0) 



(14) 



where 'ypec(O) = 300 km s~^ is the present-day galaxy veloc- 
ity dispersion (Peebles 1993). We assume that this velocity 
is directed along the line of sight. Table 4 shows the value of 
the optical depth of the three lowest rotational transitions 
of LiH as function of time, computed according to eq. (A9). 

Anisotropics of the CBR generated by the peculiar mo- 
tion of the protogalactic cloud can be positive or negative 
depending on whether the cloud is approaching of receding 
(see eq. [1]). Although the total peculiar velocity of a col- 
lapsing cloud results from the combination of Vpec and the 
collapse velocity Vcoii, we see from the values listed in Ta- 
ble 3 that the latter can be neglected for most of the collapse 
phase. For a perturbation evolving in a homologous way, 
the instant of turn around, when the radial collapse veloc- 
ity exactly compensates the overall cosmological expansion, 
presents interesting peculiarities. At this time, in fact, all the 
molecules in the cloud are "at rest" , and they all contribute 
to the optical depth. Similar considerations hold for the line 
width: at the redshift of turn-around the width of a line is the 
thinnest, practically equal to the thermal width (Ai//i^)th. 
This effect was originally predicted by Zeldovich (1978), who 
pointed out that weak and narrow spectral features gener- 
ated in protogalaxies (or protoclusters) reaching their radii 
of turn-around could be superimposed to the spectrum of 
the CBR. 



5 DISCUSSION AND CONCLUSIONS 

The presence of LiH in collapsing protogalactic clouds may 
induce fluctuations of the CBR brightness temperature on 
scales of ^ 10" which can be in principle revealed by the use 
of arrays of telescopes and aperture synthesis techniques. 
However, it is clear from the values of tliH listed in Table 4 
that the expected spectral features of LiH in protoclouds are 
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J=l-0 3=2-1 3=3-2 
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Figure 7. Relation between tlie angular size 5 of a 10^'^ Mq 
collapsing cloud and the redshifted frequency u of the three lowest 
rotational transitions of LiH. Boxes represent indicatively regions 
of the plane ii-O that can be explored with current teloscopes: 
IRAM Plateau-de-Bure Interferometer (a), IRAM 30 m (b and 
c). 



extremely weak: the associated anisotropics in the CBR tem- 
perature, according to eq. (1) are completely negligible dur- 
ing most of the collapse phase. Although the optical depth 
rises very rapidly towards the end of the collapse phase, the 
resulting value of AT/T cannot be predicted quantitatively 
on the basis of the simple model adopted in this work, and 
a more realistic approach should be taken. 

The relation between the angular size 9 of the model 
cloud and the redshifted frequencies of the three lowest 
rotational transitions of LiH is shown in graphic form in 
Fig. 7, where the coverage of the plane v-O of current tele- 
scopes is also schematically displayed. Probably the most 
appealing observational prospective is to search for CBR 
temperature fluctuations associated with the J = 1-0 rota- 
tional transition of LiH in "evolved" protogalaxies with the 
IRAM Platcau-dc-Burc Interferometer at 3 mm. Although 
the actual level of temperature fluctuations is hard to pre- 
dict quantitatively for such evolved objects, it might be close 
to the sensitivity limit of this instrument. It is also impor- 
tant to remind that the values of AT/T obtained above are 
proportional to the assumed value of the primordial abun- 
dance of Li (the results presented above are obtained with 
[Li/H]p = 2.1 X 10-1°). 

The main conclusions of our study are: 

(a) due to its low rate of radiative association, the final 
abundance of LiH is very small and the molecule is only a 
minor constitutent of the primordial gas. In our standard 
model, characterized by Ho = 67 km s-^ Mpc-^, r^io = 4.5, 
flo = 1, the abundance is /(LiH) ~ 10"" at .z = 10. Such 
a low value casts serious doubts about the use of LiH as 
a probe of the physical conditions of the gas at very high 
redshift. These findings reinforce the conclusions reached by 



Table 4. Optical depth of the first three rotational transitions of 
LiH for a 10^" Mq collapsing primordial cloud 
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SLD in a similar analysis. The abundance of LiH"*" for z ^ 18 
is about two orders of magnitude larger than the abundance 
of LiH. 

(b) Taking into account the possibility of a higher value 
of [Li]p, the final abundance of LiH and LiH^ may increase 
up to ~ 10-1^-10-1''. 

(c) The estimate of the optical depth due to resonant 
scattering of the CBR photons on LiH molecules yields a 
value of rSif ~ 10""-10-i'' at a frequency of 30 GHz. 
Most of the contribution to the optical depth comes from 
rotational transitions beteween high J levels. The attenu- 
ation of pre-existing anisotropics of the CBR is therefore 
completely negligible. 

(d) The enhancement of the LiH abundance during the 
collapse of protogalactic clouds gives rise to fluctuations in 
the temperature of the CBR associated with the peculiar 
motion of the object. For a lOi" Mq protogalactic cloud 
reaching its turn-around at 0ta = 10, the expected level of 
anisotropy is far too low to be of any observational relevance 
during most of the collapse phase. The actual abundance of 
LiH in the latest phases of collapse cannot be quantitatively 
predicted by the methods presented here. 
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APPENDIX A: RADIATIVE TRANSFER IN 
THE EXPANDING UNIVERSE 

The equation of radiative transfer in the expanding Universe 
can be written as (e.g. Peebles 1971) 



I dJ^ _ 

C dz ~ Hoil + 2)2V1 + ^02 ' c{l + z)' 

where 

= rrlllAul 1 4>W ^ '^u')^ 

]u = —nuAuicpiy - Vui)- 

Here 4>{i> — i>ui) is the normalized line-shape function for a 
transition of frequency Vui- 

In terms of the function Ji, = Ji,{l + , eq. (Al) 
takes the form 



(Al) 



(A2) 



(A3) 



CKv dz (1 + a)^ 

where Sv is the source function defined as 



g ^ '^hiy'^ I gull _ ^ 



(A4) 



(A5) 



Eq. (A4) can be formally solved in the usual way multiplying 
each term for the integrating factor where 



I Ho{l 



+ zyvi + ^oz 



dz, 



obtaining 



+ e~ 



(1 + 2)3 



'J^dr' 



(A6) 



(A7) 



Consider now a single transition of frequency Uui- Ab- 
sorption of radiation observed today at frequency v occurs 
in an interval of the order of the thermal width Ai^th around 
the Une frequency fui, or, equivalently, in a range of redshift 
Azi = (1 + 2i)Afth/^'ui around the "interaction" redshift Zi 
defined by the condition 



!/(l + Zi) = l^ul 



(A8) 



If the Doppler width of the line is small, as it is the case 
in the present context, then the size (in redshift) of the in- 
teraction region around Zi is also small, and the integral in 
eq. (A6) can be simplified as 



r(2i) = 



where 



-A ,^ 
li 91 



1 - 



giriuizi) 



guni{zi) 



ni{zi) 



(A9) 
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(AlO) 



is the "cosmological" velocity gradient, due to the expansion 
of the Universe. The region contributing to the optical depth 
has a size 



A ,1 . 
A£ = c—Azi 
dz 



Vui \dv/dr\:, 
In the same way, since 

" dT — 



(All) 



(A12) 



eq. (A7) gives 

Mr(zi)] = J.(0)e-"("') + e-"("')[l - e-"'"'>]5.(;Zi), (A13) 

and the relative perturbation in the CBR intensity at the 
present epoch is 



J. 
where 

R{zO 



= [R{zO-l][l-e-^^^'^] 



guni(zi) 
giriuizi) 









jexp 


-'} 


_kT,{zi) 



(A14) 



(A15) 



For a scries of transitions of different frequencies i^'*'-', 
the total optical depth at any observed frequency v is given 
by the sum r = r'*'' of all the single-line contributions 
t'*^ evaluated at the rcdshift of interaction zp' — v'^^^ /ii. 

Remember that for black-body radiation the relative 
distortion in the radiation intensity Jv, and in the brightness 
temperature T are related by 

AJ. xe^ AT 
where x = hv/kT. 



APPENDIX B: RADIATIVE AND 
COLLISIONAL PROCESSES 

The calculation of the level populations of LiH requires the 
knowledge of radiative and collisional rate coefficients [Rij 
and Cij) from each state i to any other state j. We have in- 
cluded in our computations only rotovibrational transitions 
satisfying the selection rules Aw = 0, ±1 for v = 0-10 and 
AJ = 0,±1 for J = 0-25. 



Bl Radiative processes 

Complete tabulations of the Einstein cocfficionts for vibro- 
rotational transitions of LiH are not available in the liter- 
ature. For the sake of completeness and homogeneity, we 
have computed all the radiative transition probabilities by 
using the appropriate formulae valid for diatomic molecules 
in the electric dipolc approximation. The rate of sponta- 
neous emission for pure rotational transitions is given by 
(e.g. Shu 1991) 

4,, 3 

(Bl) 



_ 6477Xi 2 



Shc^ " 2 J + 1 ' 

where Do = 5.888 debyes (Zemke & Stwalley 1980) is the 
value of the dipole moment evaluated at the equilibrium sep- 
aration for the two nuclei. For roto-vibrational transitions, 



A(v,J)^(v-l,J±l) = 

where 

J/(2J + 1) 



87r' 



fiuo ( dr )o^'^"'' 



if J- 
(J+l)/(2J+l) ifj- 



J-1, 
J + 1, 



(B2) 



(B3) 



1^0 = 1405.65 cm"^ (Huber & Herzberg 1979) is the fun- 
damental frequency of oscillation of LiH, and (dD/dr)o = 
0.4132e (Zemke & Stwalley 1980) is the derivative of the 
dipole moment evaluated at the equilibrium position. When 
a comparison is possible, the Einstein coefficients computed 
as indicated above agree with the accurate but sparse values 
tabulated by Zemke & Stwalley (1980) to better than 10%. 

The radiative transitions rates Rij can be expressed as 
function of the corresponding Einstein coefficients Aij 



_. A,,[1 + I3,,{T,)] iii>j 
^ {9i/9i)AiiPii{T.) ifi<i 



where 



(Ei-Ej\ 1-1 



(B4) 



(B5) 



B2 Collisional processes 

We consider only excitation and de-excitation of LiH by col- 
lision with H atoms, and we factorize dj — n(H)7ij. The 
cross sections aoj for the excitation of rotational levels of 
LiH by collision with He atoms have been calculated by Jen- 
drek & Alexander (1980) on the basis of the potential sur- 
face for the (LiH, He) system determined by Silver (1980). 
These theoretical cross sections are in good agreement with 
experimental values (Davies 1986). The collisional excita- 
tion rates 7oj(Tg) can be obtained from the known cross 
sections by integrating over a a maxwellian distribution of 
particle velocities. From the cross sections given by Jendrek 
& Alexander (1980) we have obtained accurate values of 
7oj at Tg = 3000 K. In order to compute upward rates for 
generic J and J' we have made use of the infinite order sud- 
den factorisation formula (Goldfiam et al. 1977) in the form 
given by Varshalovich & Khersonskii (1977): 



7jj'(rg) = (2j + i)^ 
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where the summation index £ assumes integer values be- 
tween (.1 — \ J' — J\ and (.2 ~ J' + J, and the term in paren- 
thesis is a Wigner 3-J symbol. This expression is valid for 
collision energies much larger than rotational energy spac- 
ings. 

The corresponding de-excitation rates arc then com- 
puted by using the principle of detailed balance. The re- 
sulting collisional rate coefficients are of the order of 1- 
2x10"^° cm^ s-i at Tg = 3000 K. Since deexcitation cross 
sections tend to be nearly constant at low energies for colli- 
sions between neutral particles (Spitzer 1978), approximate 
values of 7oj(rg) at lower temperatures can be obtained by 
multiplying the de-excitation coefficients at Tg = 3000 K by 
the factor (Tg/3000 K)^/^. An additional factor \/3 must be 
adopted to account for the difference in the reduced masses 
of the (Li,H) and (Li,He) systems. 
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